Heavy particle concentration in turbulence at dissipative and inertial scales 
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Spatial distributions of heavy particles suspended in an incompressible isotropic and homogeneous 
turbulent flow are investigated by means of high resolution direct numerical simulations. In the 
dissipative range, it is shown that particles form fractal clusters with properties independent of the 
Reynolds number. Clustering is there optimal when the particle response time is of the order of the 
Kolmogorov time scale t v . In the inertial range, the particle distribution is no longer scale-invariant. 
It is however shown that deviations from uniformity depend on a rescaled contraction rate, which 
is different from the local Stokes number given by dimensional analysis. Particle distribution is 
characterized by voids spanning all scales of the turbulent flow; their signature in the coarse-grained 
mass probability distribution is an algebraic behavior at small densities. 

PACS numbers: 47.27.-i, 47.10.-g 



Understanding the spatial distribution of finite-size 
massive impurities, such as droplets, dust or bubbles sus- 
pended in incompressible flows is a crucial issue in engi- 
neering [![, planetology Q and cloud physics 0. Such 
particles possess inertia, and generally distribute in a 
strongly inhomogeneous manner. The common under- 
standing of this long known but remarkable phenomenon 
of preferential concentrations relies on the idea that, in a 
turbulent flow, vortexes act as centrifuges ejecting parti- 
cles heavier than the fluid and entrapping lighter ones [3| . 
This picture was successfully used (see, e.g. 0,0]) to de- 
scribe the small-scale particle distribution and to show 
that it depends only on the Stokes number S v — t/t v , 
which is obtained by non-dimcnsionalizing the particle 
response time r with the characteristic time t, ; of the 
small turbulent eddies. 

In this Letter, we confirm that this mechanism is rele- 
vant to describe the particle distribution at length-scales 
which are smaller than the dissipative scale of turbulence, 
77. In particular, maximal clustering is found for Stokes 
numbers of the order of unity. However, we show that 
stationary particle concentration experiences also very 
strong fluctuations in the inertial range of turbulence. In 
analogy with small-scale clustering, it is expected that for 
r 3> 77 the relevant parameter is the local Stokes number 
<S r = t/tt, where r r is the characteristic time of eddies 
of size r [6[ . Surprisingly, we present evidences that such 
a dimensional argument does not apply to describe how 
particles organize in the inertial range of turbulence. We 
show that the way they distribute depends on a scale- 
dependent rate at which volumes are contracted. 

In very dilute suspensions, the trajectory X{t) of small 
spherical particles much heavier than the fluid evolve ac- 



cording to the Newton equation [7] 

tX = u(X,t) X 



(1) 



The response time r is proportional to the square of the 
particles size and to their density contrast with the fluid. 
Here, we neglect buoyancy and particle Brownian diffu- 
sion. As stressed in [8|, [9(, diffusion may affect concen- 
tration of particles. However we assume that the typical 
lengthscale below which particle diffusion becomes dom- 
inant is much smaller than the Kolmogorov scale rj of 
the fluid flow, and than any observation scale considered 
here. In most situations, particles are so massive that 
such an approximation is fully justified. 

The fluid velocity u satisfies the incompressible Navier- 
Stokes equation 



d t u + u ■ Vit = — Vp- 



vV 2 u 



F with V u = 0, (2) 



p being the pressure and v the viscosity; F is an ex- 
ternal homogeneous and isotropic force that injects en- 
ergy at large scales L with a rate e = (F ■ u). We per- 
formed direct numerical simulations of |(2J) by means of 
a pseudo-spectral code on a triply periodic box of size 
C = 2ir with 128 3 , 256 3 and 512 3 collocation points cor- 
responding to Reynolds numbers (at the Taylor micro- 
scale) R\ ~ 65, 105 and 185, respectively. Once the flow 
u is statistically stationary, particles with 15 different 
values of the Stokes number in the range S v € [0.16, 3.5] 
(N = 7.5 millions of particles per S v ), are seeded homoge- 
neously in space with velocity equal to that of the fluid, 
and are evolved according to H} for about two large- 
scale eddy turnover times. After this time, particle mass 
distribution reaches a statistical steady state too. Mea- 



FIG. 1: (Color online) The correlation dimension T>2 vs S v 
for the three different R\. Also shown the probability V to 
find particles in non-hyperbolic (rotating) regions of the flow, 
for Ke\ — 185 (multiplied by an arbitrary factor for plotting 
purposes). T>i has been estimated taking into account also 
subleading terms, as described in [l2l ]. 



FIG. 2: (Color online) (a) The modulus of the pressure gra- 
dient, giving the main contribution to fluid acceleration, on a 
slice 512x512x4. B/W code low and high intensity, respec- 
tively. Particle positions in the same slice are shown for (b) 
5, ; =0.16, (c) 0.80 and (d) 3.30. Note the presence of voids 
with sizes much larger than the dissipative scale. 



surements are then performed. Details on the numerics 
and on the transient are reported in ■ 

Below the Kolmogorov length-scale rj where the veloc- 
ity field is differentiable, the motion of inertial particles 
is governed by the fluid strain and the dissipative dy- 
namics leads their trajectories to converge to a dynam- 
ically evolving attractor. For any given response time 
of the particles, their mass distribution is singular and 
genericall y sc ale-invariant with fractal properties at small 
scales 
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In order to characterize particle clusters 
at these scales, we measured the correlation dimension 
I?2, which is estimated through the small-scale algebraic 
behavior of the probability to find two particles at a dis- 
tance less than a given r: P2M ~ r V2 . The dependence 
of T>2 on S v and Re\ is shown in Fig. [T] Notice that T>2 
depends very weakly on Re\ in the range of Reynolds 
numbers here explored. A similar observation was done 
in [f| , where particle clustering was equivalently charac- 
terized in terms of the radial distribution function. This 
indicates that r, ( , which varies by more than a factor 2 
between the smallest and the largest Reynolds numbers 
considered here, is the relevant time scale to characterize 
clustering below the Kolmogorov scale rj. For all values 
of i?e a maximum of clustering (minimum of T)<i) is ob- 
served for S v ~ 0.6. For values of S v larger than those 
investigated here, T>2 is expected to saturate to the space 
dimension [12| . For small values of S v , particle positions 
strongly correlate with the local structure of the fluid ve- 
locity field. A quantitative measure is also given in Fig. [I] 
by the probability V to find particles in non-hyperbolic 
regions of the flow, i.e. at those points where the strain 
matrix has two complex conjugate eigenvalues. Note that 
V attains its minimum for values of S v close to the min- 
imum of T>2- This supports the traditional view relating 
particle clustering to vortex ejection. As qualitatively ev- 
idenced from Fig. [21 particle positions correlate with low 



values of the fluid acceleration. This phenomenon was 
already evidenced in [l(| where a statistical analysis of 
the fluid acceleration at particle positions was performed, 
and also in [l3| for 2D turbulent flows. 

From Fig. [2l it is clear that fluctuations in the par- 
ticle spatial distribution extend to scales far inside the 
inertial range; this confirms the experimental measure- 
ments of [14[. Moreover, for Stokes numbers of the or- 
der of unity (Figs. [2t and[5Ji), we note that the sizes of 
voids span all spatial scales of the fluid turbulent flow, 
similarly to what observed in 2D inverse cascade turbu- 
lence dEm. To gain a quantitative insight, we consider 
the Probability Density Function (PDF) P rjT (p) of the 
particle density coarse-grained on a scale r inside the in- 
ertial range, that is the probability distribution of the 
fraction of particles in a cube of size r, normalized by 
the volume r 3 of the cube. For tracers, which are uni- 
formly distributed, and for an infinite number of parti- 
cles N — > 00, this PDF is a delta function on p = 1. 
For finite N, the probability to have a number n of par- 
ticles in a box of size r is given by the binomial distri- 
bution and leads to a closed form for the PDF of the 
coarse-grained density p = nC 3 /(Nr 3 ). In our settings, 
the typical number of particles in cells inside the inertial 
range is several thousands. Thus, finite number effects 
are not expected to affect the core of the mass distri- 
bution p — O(l), but they may be particularly severe 
when p < 1 because of the presence of voids. To reduce 
this spurious effect, we computed the quasi-Lagrangian 
(QL) mass density PDF P^,t -{p)i which corresponds to 
a Lagrangian average with respect to the natural mea- 
sure [16| , and is obtained by weighting each cell with the 
mass it contains. For statistically homogeneous distribu- 
tion, QL statistics can be related to the Eulerian ones by 
noticing that {pDql = if% +1 ) (see, e.g., 
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FIG. 3: (Color online) The quasi-Lagrangian PDF of the 
coarse-grained mass density p r in log-log scale for S v = 0.27, 
0.37, 0.58, 0.80, 1.0, 1.33, 2.03, 3.31 (from bottom to top) at 
scale r = £/16 = 0.39. The dashed line represents a uniform 
distribution. Inset: exponent of the power law left tail a vs 
S v - It has been estimated from a best fit of the cumulative 
probability which is far less noisy. Data refer to Re\ = 185. 



Figure [3] shows P r .V (p) , for various response times 
t and at a given scale (r = C/16) within the inertial 
range. Deviations from a uniform distribution can be 
clearly observed, they become stronger and stronger as 
r increases. This means that concentration fluctuations 
are important not only at dissipative scales but also in 
the inertial range. A noticeable observation is that the 
low-density tail of the PDF (which is related to voids) 
displays an algebraic behavior Pr® L \p) ~ p a ( r > T \ This 
means that the Eulerian PDF of the coarse-grained mass 
density has also a power-law tail for p <C 1, with expo- 
nent a — 1. The dependence of this exponent for fixed r 
and varying the Stokes number S v is shown in the inset. 
For low inertia (S n — > 0), it tends to infinity in order 
to recover the non-algebraic behavior of tracers. At the 
largest available Stokes numbers, we observe a — » 1, in- 
dicating a non-zero probability for totally empty areas. 
Note however that a is expected to go to infinity when 
S-q — ► oo with r fixed, since a uniform distribution is ex- 
pected for infinite inertia. Algebraic tails are relevant to 
various physical / chemical processes involving heavy par- 
ticles. The particle distribution in low density regions is 
an important effect to be accounted for, e.g., modelling 
the growth of liquid droplets by condensation of vapor 
onto aerosol particles |8[ and of aerosol scavenging [3j . 

Fixing the response time r and increasing the obser- 
vation scale r reproduces the same qualitative picture as 
fixing r and decreasing r. A uniform distribution is re- 
covered in both limits r — + oo or r — » 0. These two limits 
are actually equivalent. At length-scales r>r; within the 
inertial range, the fluid velocity field is not smooth: ac- 



cording to Kolmogorov 1941 theory (K41), velocity incre- 
ments behave as S r u ~ (er) 1 / 3 . K41 theory suggests that 
fluctuations at scale r are associated to time scales of the 
order of the turnover time t t — r/S r u ~ £ _1 / 3 r 2 / 3 . This 
implies that for any finite particle response time r and at 
sufficiently large scales r, the local inertia measured by 
S r = t /r r becomes so small that particles should behave 
as tracers and distribute uniformly in space Devia- 
tions from uniformity for finite S r are expected not to 
be scale-invariant Q . In particular, observations from 
random (^-correlated in time flows [171 ] suggest that par- 
ticle distribution should depend only on the local Stokes 
number S r . However, as explained below, this argument 
does not seem to apply to realistic flows. 

The presence of inhomogeneities in the spatial distri- 
bution of particles is due to a dynamical competition be- 
tween stretching, folding, and contraction due to their 
dissipative dynamics. At small scales where the flow 
is spatially smooth, these mechanisms are described by 
the Lyapunov exponents and their finite-time fluctua- 
tions, associated to the growth and contraction rates 
of infinitesimal volumes. At larger scales, competition 
between particle spreading and concentration is also at 
equilibrium. The spatial distribution of particle at a 
given scale r ^> r\ should only depend on the time scale 
given by the inverse of the contraction rate j r ^ T of a size- 
r blob B r of particles with response time r. However in 
the inertial range, the flow is not differentiable and an 
approach based on Lyapunov exponents cannot be used. 
In particular, the rate 7 rjT should depend on r. To es- 
timate this contraction rate in the limit of small r, we 
make use of Maxey's approximation [2l[ of the particle 
dynamics Q 

X w v(X, t) with v = u — t (d t u + u ■ Vu) , (3) 

meaning that particle trajectories can be approximated 
by those of tracers in the synthetic compressible veloc- 
ity field v. The inertial correction is proportional to the 
fluid acceleration that, in turbulent flows, is itself domi- 
nated by pressure gradient. The contraction rate of the 
blob B r with volume r 3 is j ryT — (l/r 3 )J ts d 3 xW -v(x). 
As V-v ~ tV 2 p, we have 7 rjT ~ (r /r 2 ) S r p, where S r p 
denotes the typical pressure increment at scale r. This 
means that the scale dependence of the contraction rate 
is directly related to the scaling properties of the pressure 
field. Dimensional arguments suggest that 5 r p ~ (er) 2 / 3 , 
so that the contraction rate scales as j r , T ~ -re 2 / 3 /r 4 / 3 . 
However, as stressed in [l9], K41 scaling for pressure is 
observable only when the Reynolds number is tremen- 
dously large (typically for Re\ > 600). In our simu- 
lations, where Re\ < 200, pressure scaling is actually 
dominated by random sweeping (20| . and we observe 
S r p ~ U (er) 1 / 3 (as in other simulations at comparable 
Reynolds, see [19(). This implies that the contraction 
rate is j r;T ~ re^-^U/r 5 ^ 3 . 
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FIG. 4: (Color online) Pdf of the coarse-grained mass in the 
inertial range for three values of the non-dimensional contrac- 
tion rate ■y r . T T v . From bottom to top: 7 r . T = 4.8 x 10 _4 /r^ 
(different curves refer to S n = 0.16, 0.27, 0.37, 0.48), 7 r , T = 
2.1 x W~ 3 /t v (for Sr, = 0.58, 0.69, 0.80, 0.91, 1.0), y r ',-r = 
7.9 x 1Q- 3 /t v (for S v = 1.60, 2.03, 2.67, 3.31). Inset: devi- 
ation from unity of the first-order QL moment for r within 
the inertial range. For comparison, the behavior oc 7^/ T 5 ob- 
tained when assuming uniformly distributed pointlike clusters 
of particles is shown as a solid line. Data refer to Re\ = 185. 

Figure |4] shows Pr® L \p) for three choices of 7 r)T ob- 
tained from various sets of values of r and r. The collapse 
of the different curves strongly supports that the distri- 
bution of the coarse-grained mass density depends only 
upon the scale-dependent contraction rate 7 r ,r- In partic- 
ular, as represented in the inset of Fig. [4j the deviations 
from unity of the first moment of the distribution col- 
lapse for all S v investigated and all scales inside the iner- 
tial range of our simulation. This quantity is the same as 
the Eulerian 2nd-order moment, giving the probability to 
have two particles within a distance r. We note that par- 
ticle distribution recovers uniformity at large scales very 
slowly, much slower than if particle were distributed as 
Poisson point-like clusters for which (p 2 ) — locr~ 3 oc 7, 9 i (- 5 
(also shown in the inset for comparison). 

We presented numerical evidence that at moderate 
Reynolds number the distribution of heavy particles at 
lengthscales within the inertial range is fully described 
in terms of a scale-dependent volume contraction rate 
7 r T oc r/r 5 / 3 . However we expect that j r ^ r oc r/r 4 / 3 at 
sufficiently large Reynolds numbers (Re\ > 600). State- 
of-the-art numerical simulations can not currently attain 
these turbulent regimes, so that only experimental mea- 
surements of particle distribution can help to confirm 
such a prediction. We have seen that particle dynamics in 
the inertial range can be directly related to the structure 
of the pressure field (and thus of acceleration). Charac- 
terizing the distribution of acceleration is thus crucial to 
understand particle clusters, and conversely, particle dis- 
tribution could be used as an experimental tool to probe 



the spatial structure of turbulent flows. 
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